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Abstract 

We describe families of potentials which act as approximate cloaks for matter waves, i.e., for 
solutions of the time-independent Schrodinger equation at energy E, with applications to the 
design of ion traps. These are derived from perfect cloaks for the conductivity and Helmholtz 
equations, by a procedure we refer to as isotropic transformation optics. If W is a potential which 
is surrounded by a sequence {^f'}^L 1 of approximate cloaks, then for generic E, asymptotically 
in n (i) W is both undetectable and unaltered by matter waves originating externally to the cloak; 
and (ii) the combined potential W + V® does not perturb waves outside the cloak. On the other 
hand, for E near a discrete set of energies, cloaking per se fails and the approximate cloaks support 
wave functions concentrated, or almost trapped, inside the cloaked region and negligible outside. 
Applications include ion traps, almost invisible to matter waves or customizable to support almost 
trapped states of arbitrary multiplicity. Possible uses include simulation of abstract quantum 
systems, magnetically tunable quantum beam switches, and illusions of singular magnetic fields . 
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Introduction. A fundamental problem is to describe the scattering of waves at energy E 
by a potential V(x), as governed by the time-independent Schrodinger equation; the related 
inverse problem consists of unique determination of V from scattering data or from boundary 
measurements. In this Letter, for a generic energy E, we construct families {V^}^ =1 of 
bounded potentials which act as an approximate quantum cloaks: for any potential W whose 
support is surrounded by the support of the scattering amplitude of W + goes to 
zero as n — > oo, so that W is, asymptotically in n, undetectable by waves at energy E. 
The central potentials are supported in a spherical annulus {t\ < r < r 2 } in M. 3 and 
have spatial oscillations of increasing amplitude as r \ r\ and decreasing quasiperiod as 
n — > oo, so that their local (resp., long range) effect on wave propagation becomes stronger 
(resp., weaker) as n — > oo. For generic E, the potential W + V® has a negligible effect on 
matter waves originating outside of its support. Alternatively, for E close to special values, 
allows the core {r < r\\ to support almost trapped states and be used to form traps for 
ions (here denoting any charged particles), almost invisible to external matter waves. An 
approximate version of the dichotomy regarding ideal cloaks for the Helmholtz equation jH 
Thm. 1] holds: If E is sufficiently separated from all eigenvalues, then with high probability 
the approximate cloak keeps particles of energy E from entering the cloaked region, see 
Figs. l(r,red) and 2(1); on the other hand, for E close to an eigenvalue, the cloaked region 
supports almost trapped states, accepting and binding such particles, leading to a new type 
of ion trap, cf. Figs. l(r,blue) and 2(r) and formula Q. 

Recently, Zhang, et al., [Ij, using ideas from transformation optics [21 E], described an 
ideal quantum mechanical cloak at any fixed energy E and proposed a physical implementa- 
tion. The construction starts with a homogeneous, isotropic mass tensor mo and potential 
Vq = 0, and subjects this pair to the singular transformation ([I]) below. The resulting m, V 
yield a Schrodinger equation that is in fact the Helmholtz equation for the corresponding 
singular Riemannian metric and thus covered by the analysis of cloaking for the Helmholtz 
equation in jlj. The treatment there shows that the potential within the cloaked region is 
undetectable by exterior measurements, and that the perfectly reflective Neumann bound- 
ary condition automatically holds for finite energy distributional solutions at the inside of 
the cloaking surface (the inner surface of the cloak), S. The cloaking mass tensor m and 
potential are both singular, and m infinitely anisotropic, at S. These combine to make 
such a quantum cloak challenging to construct, with ultracold atoms in an optical lattice 
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proposed in [T] possible realization. 

The approach in this Letter is to forgo the perfect functioning of the ideal quantum cloak, 
and to describe bounded potentials which function as approximate cloaks with respect to 
the isotropic mo, with the failure to cloak perfectly in fact a controllable feature that can be 
taken advantage of for applications described below. The are found by means of isotropic 
transformation optics, a technique we introduce for avoiding the singular and anisotropic 
behavior, difficult to physically realize, of material parameters that commonly occur in 
transformation optics-based designs [SI El E] ; more details and proofs can be found in [7] . 

Inverse scattering and conductivity cloaks. There is an enormous literature on unique de- 
termination of a potential V from scattering of plane waves at energy E by the Schrodinger 
equation (—V 2 + V)ip = Eip, as encoded in the scattering amplitude ay(E, 9, uj). Equiv- 
alently, for compactly supported V one may consider the near-field measurements of wave 
functions at the boundary of a large region Q, as encoded in the Dirichlet-to- Neumann (DN) 
operator, Av(E)(i(}\qq) = d u ip\ssi where v is the normal of dfl [8]. Ideal quantum cloaking 
can be considered as giving highly singular examples of nonuniqueness, but in order to 
construct approximately cloaking potentials V^, we first need to recall the ideal cloaking 
conductivity a\ of [9]. For simplicity, we describe the cloak on B3 — B1, with Br denoting the 
central ball of radius R in M 3 , so that the cloaking surface, the interface between the cloaked 
and uncloaked regions, is £ = {r = 1}. Subjecting a homogeneous, isotropic conductivity 
cr to a singular transformation ("blowing up a point"), we constructed [S] certain singular, 
anisotropic conductivity tensor fields on B 3 — B± which, when augmented by any nonsingular 
(bounded above and below) conductivity on B±, result in conductivities on B% giving the 
same electrostatic boundary measurements as <To. (For related results in dimension two, see 
[TOf ITT].) The same construction, applied to the electric permittivity and magnetic perme- 
ability rather than the conductivity tensor, was used [3] to propose an electromagnetic (EM) 
cloak at nonzero frequency, and a microwave experimental realization of a variant of this 
design reported in (T3j. Ray-based cloaking for dimension two was proposed in [5], while 
potentials transparent for rays are in [12J. 



To be precise, let F = (F 1 , F 2 , F 3 ) : B3 — {0} — > B% — B\ be the singular transformation, 
for r e M 3 , 



which results in the transformed version of ctq on B 3 — Bi, augmented for simplicity by 2o"o 




(1) 
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on Bi, 

a 1 = F,a , r e B 3 - B 1 ; a x = 2 a , r e B v (2) 

Here, denotes the change-of- variables action of F on tensors a = (<J^ h ), writing r = 
(x 1 , x 2 , x 3 ), 

[ * ' [ ' detm\ ^ dxv dx^ 

ldx k J p,q=l T=F -itt 

This a i has a singularity at E, both in that one of the eigenvalues (corresponding to the 
radial direction) tends to as r \ 1, and that there is a jump discontinuity across S, 
within which a\ is non-singular. Aside from the radius of the outer ball and the factor 
2 in the second formula of ([2]), o\ is the conductivity introduced in [9] and shown to be 
indistinguishable from <7 > vis-a-vis boundary measurements at dB 3 of electrostatic fields. 
Consider also the corresponding acoustic equation, 

d, (a[ j d jU ) + Ea lU = 0, a x = (det^])" 1 . (3) 



Then, using boundary measurements at dB% of the acoustic waves of frequency yE, the 
pair oi, ai is indistinguishable, from a , a = 1 jU dH [T5J US]. 

In the case of ideal cloaking, the waves u within the cloaked region have a simple descrip- 
tion [U Thm.l]: (I) if E is not a Neumann eigenvalue of — V 2 + W in the cloaked region, then 
u must vanish there; or (II) if E is an eigenvalue, then u can be an associated eigenfunction 
there, while possibly vanishing outside of the cloaking surface. The Dirichlet eigenvalues 
and eigenfunctions of |3) on B 3 can be separated into (Ep 1 ,^) and {E~' 1 , uj) with and 
uj supported in B 3 — B\ and Bi, resp. 

Approximate cloaks and failure of cloaking. For 1 < R < 2, let or be given by the 
same formulae as in (|2]), but applied on B 3 — Br and Br, resp., and similarly for or in (|3]). 
Observe that for each R > 1, an and a# are nonsingular, however, their lower (resp., upper) 
bounds go to (resp., oo) as R \ 1. (Similar truncations of EM cloaks have been studied 

m puna eg.) 

When o\ y a\ are replaced by <tr, or there is a decomposition similar to the one above for 
(ij), with eigenvalues and eigenfunctions (Ep R , vf) and (EJ ,R , vj) concentrating in B 3 — B\ 
and Bi, resp., and Ef' R converges to E-' 1 as R — » 1. The solution of the boundary value 
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problem <9, (o~ R djv) + Ea R v = on B 3 , v\qb 3 = f, has an eigenfunction expansion 



An approximate version of dichotomy (I)-(II) holds for approximate cloaks: When E is 
not equal to any Ep 1 , one can show the DN operators for the <7r, Or converge to that for 
<Ji,ai as R — > 1; physically, this means that the boundary measurements of pressure and 
the normal component of the particle velocity for the approximate cloaks tend to those for 
the ideal cloak, which are themselves the same as for do, oq. However, if E is close to some 
Ep R , the corresponding term in (4J) may dominate the others, in which case the solution 
v, having a large coefficient of vj, concentrates in B\. Since v~ cannot vanish identically 
in B 3 — Bi, both the near- field measurements on the boundary dB 3 and far-field patterns 
differ from those corresponding to (Tq, a . This interior resonance corresponds to a sound 
wave almost trapped within the cloak. 

Isotropic transformation optics. A well known phenomenon in effective medium theory 
is that homogenization of isotropic material parameters may lead, in the small-scale limit, 
to anisotropic ones [19]. We exploit this, using ideas from [2Ql E3], to approximate the 
anisotropic, almost cloaking o~r by isotropic conductivities a Rte so that for e > the o~R t6 , a R 
also function as approximate acoustic cloaks [7]. The cj^ e (r) are spatially highly oscillating, 
with values varying roughly between the extreme eigenvalues of or near r. (In the context 
of EM cloaking, thin concentric layers of homogeneous, isotropic media were considered in 
[221 [23].) 

Approximate Schrddinger cloaks. The gauge transformation if) = \fou reduces the 
acoustic equation (|3|, with nonsingular isotropic conductivity a = a^ e in place of the 
anisotropic o\ and or in place of a%, to the Schrddinger equation at the same energy E, 
(- V 2 + Vl) if) = Eif>, where V£ = V 2 ^ e )/^, e + E (l - afa^) ■ As a R>e is highly 
oscillatory, V Re consists of a quasiperiodic pattern of concentric radial potential barriers and 
wells of increasing amplitudes and decreasing widths as e \ 0. The radial profile of the 
potential over one quasiperiod (i.e. in one spherical layer) is in Fig. 1 (left) . Since <7r,<e = 1 
near dB 3 , the boundary measurements of solutions of these Schrddinger equations at dB 3 
coincide with those for the corresponding acoustic equations. By the convergence for the 
acoustic equations, we can choose R \ 1, e \ 0, so that the boundary measurements for 
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these Schrodinger equations converge to those for the acoustic equation ([3]) at energy E, 
which are the same as for the Schrodinger equation in free space. Our main result is: 

Approximate Quantum Cloaking. Let W be a bounded potential on B\, and E 6R be 

neither a Dirichlet eigenvalue of the free Hamiltonian — V 2 on B% nor a Neumann eigenvalue 
of — V 2 + W on B\. Then there exists a sequence of cloaking potentials on B 3 such that 
the DN operators A w+v e(E) — > A (E) as n — > oo. I.e., as n — > oo the potential W + V® is 
indistinguishable at energy E from the zero potential by near-field measurements; a similar 
result holds for far-field patterns. W is thus approximately cloaked when surrounded by ■ 
As any specific measurement device has a limited precision, this means that it is possible 
to design a potential to cloak an object within from any single-particle measurements made 
using that device at energy E. 

Numerics: We use analytic expressions to compute the wave function ip for an incident 
plane wave with ipi nc (x) = ae lkv ' d . The computations are made without reference to physical 
units, using a = 1,E = 0.5, k = yE. The cloak is based on R = 1.005, corresponding 
to an anisotropy ratio of cr^ at S/j = {r = R} of 4 x 10 4 . In the simulations we use a 
cloak consisting of 20 homogenized layers inside and 30 homogenized layers outside of the 
cloaking surface Inside the cloak we have located a centrally symmetric step potential, 
W(x) = CinnX[o,o.9]{ r ) ■ in the simulations, the cloaking potential V® and the energy E are 
the same in all figures, but we vary the constant Cj nn . In Fig. l(r, red) and Fig. 2(1) we have 
Cinn — —98.5, and ip is the wave in IR 3 produced by an incoming plane wave. In Fig. l(r, 
blue), with c inn = +1.858, and in Fig. 2(r), with c inn = —71.45, there is no incoming wave, 
but an almost trapped state corresponding to an excited state in the cloaked region. 

Applications: Almost trapped states and ion traps. A version of the dichotomy for 
approximate acoustic cloaks described above also holds for approximate quantum cloaks, 
since the solutions differ only by a gauge transformation. As a consequence, given an energy 
E, the approximate cloak may be such that E either is or is not an eigenvalue for the ideal 
model. This results in B\ becoming either (J) an almost cloaked region that with a high 
probability does not accept energy E particles from outside S; or (II) a trap that supports 
almost trapped states, which correspond to a particle at energy E trapped in B\ with high 
probability. A design of either type could possibly be implemented by an array of ac and 
dc electrodes with total effective potential approximating V® for large n [21] . This leads to 



6 




r 



FIG. 1: Left: The radial profile of the potential over a typical quasi-period 1.5 < r < 1.51. 
The potential in R = 1.005 < r < 2 is obtained by repeating similar profiles, with increasing 
amplitudes as r \ R. Right: Re^ on a segment {(x, 0,0) : < x < 3} for the same cloaking 
potential and two different cloaked W's. For the red curve, -E is not close to an interior 
eigenvalue and ip is produced by an incoming plane wave. For the blue curve, E is a Dirichlet 
eigenvalue on B3 and ip is an almost trapped state. 

a new type of trap for ions, differing from, e.g., the Paul [2S]> Penning or Zajfman [27] 
traps, and justified on the level of quantum mechanics. Furthermore, the trap may be made 
tunable by including a dc electrode in the trapped region, corresponding to a Coulomb W; 
varying the charge changes whether or not E is as eigenvalue of —V 2 + W and thus which 
of (7) or (II) holds. 
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FIG. 2: Left: Heip in B% for ip resulting from an incident plane wave. E is not near an interior 
Neumann eigenvalue; the matter wave passes unaltered. Right: An excited almost trapped state. 
E = E- ' is an energy close to a Neumann eigenvalue of B\ , for which the ideal cloak supports an 
trapped state. 

Topological ion traps. The basic construction outlined above can be modified to make the 
wave function on B\ behave as though it were confined to a compact, boundaryless three- 
dimensional manifold, topologically but not metrically the three-dimensional sphere, § 3 . By 
suitable choice of metric, the energy level E can have arbitrary multiplicity for the interior 
of the resulting trap, allowing one to implement physical systems mimicking matter waves 
on abstract spaces. As the starting point one uses not the original cloaking conductivity 
o\ (the single coating construction), but rather a double coating [U Sec. 2], which we denote 
here by . This is singular from both sides of E, and in the EM cloak setting corresponds 
to coating both sides of £ with metamaterials. See [7J Fig. 7]. By [U Sec. 3.3], the finite 
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energy solutions of the resulting Helmholtz equation on B s split into direct sums of waves 
on B 3 — Bi, as for <7i, and waves on Bi which are identifiable with eigenfunctions of the 



standard metric on S 3 , then nonground states are degenerate and of high multiplicity, while 
a generic choice of g yields simple energy levels [28]. By suitable choice of g, any finite 
number of energy levels and multiplicities can be specified [29j, allowing traps supporting 
almost trapped states at energy E of arbitrary degeneracy. 

Magnetically tunable quantum beam switch. Consider a beam of ions of energy E, leav- 
ing an oven and traveling through a tube T = {0 < p < po, < 9 < 2ir, < z < L} 
(in cylindrical coordinates). Treating the ions as matter waves, place in T several almost 
trapping traps of the type described above, centered at points on the z-axis, forming a po- 
tential V = Ylf=i Vni z ~ z j)- The techniques above may be applied to magnetic Schrodinger 
equations and DN operators on a region Q, 



We design the traps so that, in the absence of a magnetic field, or for small field strengths, 
the traps act as cloaks and thus the ions pass through T unhindered. However, if a homo- 
geneous magnetic field is then applied to the tube, chosen so that the magnetic Schrodinger 
operator has E as a Neumann eigenvalue inside each trap, then there is a large probability 
that an ion passing the jth trap will bind to that trap. If N is large enough, then the 
probability that any ion traveling the length of T will become bound is ~ 1, and T thus 
functions as a magnetically controlled switch for the beam of ions. 

Magnified magnetic fields. For a homogeneous magnetic field with linear magnetic poten- 
tial A, one can obtain a sequence of electrostatic potentials V n for which lim^oo Av„,a(E) = 
A ^(E), with A singular at a point. I.e., in the presence of a homogeneous magnetic field, 
the V n produce scattering or boundary measurements that tend, as n — > oo, to those of 
the zero electrostatic potential in the presence of a magnetic field blowing up at a point, 
|A(r)| ~ r , |B(r)| ~ r~ 2 , giving the illusion of a locally singular magnetic field [7J. 
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Laplace-Beltrami operator — V 2 g 



on (§ 3 ,g) with eigenvalue y/E. If one takes g to be the 




A VA (E)(f) = d u ij\ dn + i (A-u)f ondQ. 
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